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^ I We analyze theoretically the propagation of surface plasmon polaritons 

cn I about a metallic corner with a finite bend radius, using a one-dimensional 

model analogous to the scattering from a finite-depth potential well. We 
\0 • obtain expressions for the energy reflection and transmission coefficients in 

2 . the short wavelength limit, as well as an upper bound for the transmittance. 

O '. In certain cases we find that propagation on non-planar interfaces may result 

c/3 I in lower losses than on fiat surfaces, contrary to expectation. In addition, 

. y ■ we also find that the maximum transmittance depends non- monotonously on 

j;^' the bend radius, allowing increased transmission with decreasing radius. 

^' 

Fri Structured materials which allow nanoscale control of light are necessary for achieving 

compact integrated photonic devices. While the size of standard optical components and 
^ ■ beams is typically set by the diffraction limit, low dimensional excitations such as surface- 

plasmon polaritons may be confined to dimensions much smaller than the wavelength of 
light. Surface-plasmon polaritons (SPPs), coupled modes of plasmons and photons, are 
excited when visible electromagnetic (EM) radiation couples into surface guided modes 
at metal-dielectric interfaces [1, 2]. When propagating along fiat interfaces, these are 
essentially two-dimensional (2D) waves, with an EM field intensity which peaks at the 
interface and decays exponentially into the two adjoining media. 

Recently, SPP waveguiding and bending in nano-patterned metallic films were studied 
[3]. Alternately, it was shown that EM energy may be efficiently transported by near 



^ 



field coupling in plasmon waveguides comprised of ordered arrays of metal nanopar- 
ticles [1]. Optical elements such as linear waveguides [5], mirrors, beamsplitters and 
interferometers [6] were recently demonstrated. 

Interestingly, while significant progress has been made in understanding SPP propaga- 
tion in nano-structures, certain fundamental issues pertaining to their guiding on smooth 
metallic films remain unknown. In particular, quantifying guiding and energy losses in 
SPPs propagating around bends in metal-dielectric interfaces is of great importance, as 
it should set a limit on feature size in certain plasmonic-circuit devices. Previously, the 
problems of refraction [7] and reflection of SPPs [8] at interfaces have been addressed in 
this context. In this Letter we present a study of the efficiency of SPP propagation at a 
curved metal-dielectric interface in the short wavelength limit. 
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Figure 1: Cross-sectional (a) and top view (b) of the SPP intensity; (b) also shows 
a drawing of the geometry. A metallic corner characterized by a dielectric 
constant e^ and Re[ej] < has a bend angle 6 > 90° and a finite bend radius 
R. The bend is confined to the region of space shown, with the center of 
curvature at the origin. The rest of space is occupied by a dielectric with eo- 
Axes Xi and X2 extend along the boundaries between Regions I and II and 
Regions II and III, respectively. In the a;i-a;2 plane. Regions I and III are semi- 
infinite. The system is also infinite in extent along the entire z axis. In (a) 
we illustrate the single- mode approximation developed in the text: the field 
profiles in Regions I (solid line) and II (dashed), calculated for ojR/c = 800, 
are well matched. In (b), the intensity is overlayed in grayscale, showing the 
overlap with the metal (dielectric constant ej) and the outer dielectric (eo). 
Arrows indicate incident and refiected fields in Region I, and transmitted field 
in Region III. 



The geometry we study is that of propagation about a rounded edge, as shown in Fig. 
1. SPPs are incident along the interface from Region I, and propagate counterclockwise 
through the bend in Region II, in the direction of Region III. We calculate the energy- 
refiection and transmission coefficients, as well as bend-induced radiation losses. For 



lossy metals absorption losses are also evaluated. 

Our approach exploits known expressions for the SPP fields in each region and matches 
them at the two ends of the bend. The procedure differs from related numerical tech- 
niques [9] in that we consider the SPP itself as the incident wave, with the goal of 
manipulating it as a well-defined quasi-particle in a non-trivial geometry. Favorable 
conditions for this will be seen to emerge in the short- wavelength limit, and we use this 
to arrive at analytic expressions. 

The solutions for SPP propagation on an infinite fiat surface and on cylindrical surfaces 
are known analytically [2]. On a fiat interface, SPPs at frequency uj are two dimensional 
waves, decaying exponentially into the two adjoining media, with decay constants 7^ = 
— cjejA/— l/(ej + eo)/c in the metal and 70 = u}eo\/ —l/{ei + eo)/c in the dielectric. In 
the limit Re[7i]-R ^ 1 the interference of SPPs in Regions I and III is negligible, allowing 
us to use the infinite fiat surface 2D solution in these regions. 

We construct the solution in Region II using the known solutions for SPPs propagating 
around the perimeter of an infinitely long cylindrical metal rod of radius R [2]. Here, 
the magnetic field is given by 
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where fcj = uj^/el/c and J„ is the Bessel function. The set {n} is determined by the 
metal boundary matching equation, 

1 J'ShR) _ 1 HJ^'^'iKR) 
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where ko = u;y^/c, Hn is the Hankel function of the first kind, and the prime denotes 
differentiation with respect to the argument. Assuming u real, one finds n to be complex, 
as a consequence of radiation loss and absorption in the bend. Since the wave depends 
on as exp[±m0] = exp[±(iRe[n] — lm[n])0], where is measured from the Xi axis, 
only solutions with Im[n] > are physical for damped propagation. 

Solving exactly for the transmission and refiection coefficients requires matching an 
infinite number of solutions at the boundaries along the Xi- and X2-axes separately. 
However, it is possible to render this problem tractable by a few simple approximations. 
Noting that the incident SPP carries momentum proportional to k = uj^/ eoeijieo + ej)/c, 
in the short wavelength limit uoR/c ^ 1 its angular momentum with respect to the origin 
is approximately equal to Re [/;;]/?. In Region II the solution has angular momentum 
proportional to the various n-values. Conservation of angular momentum dictates that 
the incident SPP couple predominantly to that cylindrical mode with n closest in value to 
kR. Therefore it is necessary to consider only a single term of the expansion. Formally, 
this is shown by examining the set {n} and noting that it contains an element m which 
minimizes the mismatch between the field profiles perpendicular to the surface. The 
role of angular momentum conservation in this matching problem is analogous to that 
of tangential momentum conservation in refraction at a dielectric interface. We call the 



clockwise and counterclockwise modes corresponding to m the fundamental modes. In 
the short wavelength limit of a fundamental mode n = m k, kR, and the decay rate is 
\J'm^ jB? — kf ^ 7j. Thus, Jmihr) ~ exp[7jr] near the interface, identical to the fields 
in the metal in Regions I and III. Similarly, in the dielectric Hm (kor) ~ exp[— 7or] [10]. 

The modes n ^ m have decay rates not as close to Re[7j] as that of the fundamental 
mode's. For this reason, it is possible to assume that in the short wavelength limit 
the incident SPPs couple predominantly to the fundamental modes and ignore other 
mode coupling. In order to satisfy the standard Maxwell boundary conditions it is 
therefore necessary to match only a small number of solutions at a single point on each 
axis, at a distance R from the origin. The boundary conditions are thus also satisfied 
approximately over the entire extent of the axes. As can be seen from Fig. 1, the mode 
mismatch at the boundaries may be very slight. The problem has now essentially become 
one dimensional (ID), analogous to scattering from a ID finite potential well [11]. Since 
the allowed ?7i- values are always complex, bound-state solutions in the well do not exist. 
This distinguishes the SPP on a bent surface from waveguide bends enclosed on all sides 
by infinite potential walls [12]. 

Applying the appropriate boundary conditions to the fields at the two boundaries 
results in the familiar expression for the transmittance 
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When the losses in the metal are accounted for, Im[??7,] increases with /2, such that when 
Im[?Ti]^ ^ 1 the transmittance becomes 

T ^ I6^1^^^e-2I-He. (3) 

The reflectance R may be obtained in a similar manner. In the limit uoR/c -^ oo these 
expressions become exact. 

For lossless metals, the bend-induced radiation losses are simply given by A = 1— T— R. 
Accounting for absorption in the metal we find that A now includes both radiation and 
absorption losses. We extract the radiation losses by integrating the Poynting vector S 
for unit incident flux in Region II at r ^ oo: 

e+<t>o 

S ■ rrc/0. (4) 

Since the radiation carries angular momentum with respect to the origin, the energy 
radiated into the far-field from the surface at = propagates at an angle 0o with 
respect to the xi axis, setting the lower integration limit in (4). In the short-wavelength 
limit only the amplitude of the forward-propagating mode is significant, therefore the 
radiation losses are well approximated by integrating only this mode. To obtain 0o we use 
a stationary phase approximation. The position-dependent phase is $ = koV + Re[?Ti]0, 
and the vector normal to a surface of constant phase is v(r) = V$ = koi + Re[m]/re0. 
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Figure 2: The upper bound for the transmittance, T„, plotted for a silver-air interface 
with 6 = 90°, as a function of bend radius R for wavelengths of A = SOOnm 
(dashed-dotted), A = 600nm (dashed) and A = 700nm (solid). Inset: T„ in 
grayscale as a function of R and A. 



The change in angle as the wave propagates a radial distance 6r is 50 = Kelm] / {kor'^)6r , 
giving 00 = J^ 7ie[m]/ {kor'^)dr = Re[m]/{koR)- 

We have carried out calculations for typical values of silver (e^ = —15 + zO.5) in air 
(eo = 1) when ujR/c = 800 and 6 = 90°. Ignoring the losses in the metal we find 
T = 0.997, R = 1.19 X 10"^ and P ^ 0.003. When the losses are included the results 
change drastically to T = 0.0516, R = 1.18 x 10"^ and P ^ 0.00282, indicating that most 
of the energy is lost to absorption in the metal. Comparing the latter overall absorption 
and radiation losses to the energy absorbed when SPPs propagate the equivalent arc 
distance on a flat surface, we find that propagation on a non-planar interface may result 
in lower losses. We explain this counterintuitive result using an analogous picture of 
semi-classical motion under an effective potential in a central potential field. In the 
short wavelength and large angular momentum limit the SPP fields propagating on the 
curved interface sample less of the metal volume than that available when propagating 
on a flat interface, hence the reduced absorption. 

We evaluate the accuracy of our result by examining the coupling efficiency A of a 
single mode on a flat interface to a fundamental mode m. We define this by 
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where rj 



/«"^"° |exp[-7o(r-/?)]rrfr 
0(1). We find that the condition A ^ 1 constitutes a stricter criterion 



for the validity of our approximation. When the latter holds, the incident SPP couples 
predominantly to the fundamental modes, making the approach described above self- 
consistent. For example, when 77 = 3, e^ = —15, and eo = 1, A^ = 0.002 for ojR/c = 800, 
rendering our result applicable. On the other hand, for ujR/c = 100 we obtain A^ = 0.3, 
signifying that the expression is not reliable because the coupling to modes n other than 
the fundamental can no longer be neglected. In this regime a more physical quantity is 
the upper bound for the transmittance, given from (3) by 

T„ = exp(-2Im[m]^). (6) 

Here we neglect reflections at the boundaries between the different regions, thus exclud- 
ing interference with the counter-propagating mode in Region 11. To understand why 
this is an upper bound, recall that in the wavelength range of interest, where the metal 
is not very lossy, Im[n] > Im[m]. Since the wave depends on n as exp[±m0], modes 
with large Im[n] decay rapidly. Thus, the transmission in the presence of coupling to 
non- fundamental modes does not exceed T„, and the latter is a true upper bound. Fig. 2 
is a plot of T„. A peak is clearly visible, moving to higher values of R as the wavelength 
increases. To the right of the peak, at large radii of curvature absorption losses in the 
metal dominate, and the maximum transmittance decreases with increasing radius. To 
the left of the peak radiation due to the high curvature is the dominant loss mechanism, 
leading to a rapid drop in T^- At very high curvature {R < lO/im) there is a change 
in trend, and T^ starts to increase with decreasing R (see Inset.) When calculating the 
radiation loss per arclength, we find that for this range of radii it increases slower than 
elsewhere, allowing T„ to increase even as R attains very small values. This anomalous 
behavior can be observed for all wavelengths, and is independent of the dispersion in 
the metal. 

The formalism developed here may also be used to analyze the complementary reversed 
geometry, where the metal occupies three quadrants in space, and the SPPs propagate 
around a dielectric void in it. Surprisingly, in this case we find that in the single mode 
approximation SPP propagation around the bend is non-radiative. Separate work will 
address radiation and absorption processes in this system in greater detail. 

In summary, we have analyzed the scattering of SPPs at a curved metal-dielectric 
interface in the short wavelength limit. Utilizing an analogy to a quantum mechanical 
ID finite square well we obtained the energy transmission and reflection coefficients. 
Interestingly, propagation on a curved interface may result in lower losses than at a flat 
metallic surface, due to the unique field distributions which arise in our system. An 
expression for an upper bound on the transmittance was also obtained, showing that at 
high curvature radiation is the main loss mechanism, while at low curvature material 
losses dominate. An unexpected behavior where the maximum transmittance increases 
with curvature was also observed. We explain this as an interplay between various loss 
rates in the system. These results shed new light on the mesoscopic behavior of SPPs, 
and should play an important role in the design and optimization of SPP devices. Future 
work will address SPP propagation in waveguides, splitters and interferometers. 
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